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Even in the absence of external influences the operability of a quantum computer 
(QC) is not guaranteed because of the effects of residual one- and two-body im- 
perfections. Here we investigate how these internal flaws affect the performance of 
a quantum controlled-NOT (CNOT) gate in an isolated flawed QC. First we find 
that the performance of the CNOT gate is considerably better when the two-body 
imperfections are strong. Secondly, we find that the largest source of error is due to 
a coherent shift rather than decoherence or dissipation. Our results suggest that the 
problem of internal imperfections should be given much more attention in designing 
scalable QC architectures. 
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1 Introduction 

Quantum computers (QCs) are subject to internal as well as external errors. 
External errors arising from system-environment interactions (i.e. decoher- 
ence and dissipation [1]) are perceived to be the primary obstacle to practical 
implementation of quantum computation [2]. However, even in the absence of 
a macroscopic environment the efficient operability of a QC is not guaran- 
teed. Unavoidable one- and two-body static internal imperfections can also 
be potentially harmful error sources. 
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Recently, the problem of internal errors has attracted some attention. Studies 
on the statistical properties of eigenspectra of flawed QCs showed that suffi- 
ciently strong residual two-qubit interactions cause the emergence of chaos [3, 4] 
and a consequent thermalization of the QC core [5] which is believed to be a 
serious obstacle for quantum computation. Destructive effects of chaos on QC 
operation include an incoherent mixing of the ideal (i.e. separable) compu- 
tational states of the qubits[3]. This mixing is expected to worsen with time 
resulting in an effective QC meltdown[4]. However, dynamical calculations 
in which the effects of internal flaws were modeled by random kicks showed 
better fidelity [6, 7] decay for chaotic than for integrable perturbations. These 
seemingly contradictory conclusions, i.e. destruction and stabilization of QCs 
with chaos, suggest that determining the effects of internal flaws is a problem 
of considerable complexity. 

The pathology of internal errors includes internal decoherence and dissipation 
and even coherent shifting of the QC operation. The relative importance of 
these effects and their dependence on residual internal couplings can only 
be determined via fully realistic dynamical simulations of QC operation. The 
objective of this Letter is to directly examine flawed QC operation in a specific 
architecture, to investigate the different types of internal errors that emerge, 
and to observe how these errors change with varying magnitudes of internal 
flaws. 

A suggested remedy for the errors caused by these flaws is to operate the QC 
below the chaos border. This could be achieved by employing weaker coupling. 
This can be a problematic solution[4], however. First of all, the duration of a 
gate operation is inversely proportional to gate coupling strength, and longer 
gate operations are more likely to be decohered due to the external influences. 
Secondly, even below the chaos border (i.e. in the integrable regime), the resid- 
ual two-qubit interactions are still strong enough to mix ideal logic states of 
qubits[8], and thus residual entanglement among qubits is unavoidable. In ad- 
dition there is evidence that the irreversible effects of internal decoherence are 
more harmful in the integrable regime [9, 10], although chaos is indeed harmful 
at high temperatures[ll]. 

We focus our study on the dynamics of a two-qubit subsystem located at the 
center of a 3x4 qubit 2D grid circuit. The qubits are modeled and parametrized 
based on a Josephson charge qubit QC[12] proposal. This scalable architecture 
has an extremely long external decoherence time of ~ 10~ 4 s, which in principle 
allows around 10 6 single qubit operations. However, internal imperfections for 
this design could be a limiting factor. Our two-qubit subsystem (i.e. the active 
part) performs a quantum controlled-NOT (CNOT) gate while the rest of the 
qubits are idle. We simulate the effects of flaws by adding residual one- and 
two-body imperfections in the idle part. The ten qubits neigboring the active 
ones are chosen to be idle for computational simplicity. 
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We monitor the dynamics of the CNOT gate by two traditional error quan- 
tifiers: purity and fidelity. Purity, while insensitive to unitary errors, gives a 
good measure of non-unitary errors. Fidelity, on the other hand, can detect 
both non-unitary and unitary errors. We find that non-unitary errors due 
to internal decoherence and dissipation are non-negligible even for relatively 
small numbers of bath qubits. However, we find that strong chaos in the bath 
degrees of freedom leads to complete suppression of non-unitary errors. Sur- 
prisingly, the system's fidelity is much worse than its purity would suggest. 
The unexpectedly large differences between the purity and fidelity show that 
unitary errors generated by coherent shifting are dominant in the system dy- 
namics. These large unitary errors persist even in the chaotic regime, and 
severely endanger the viability of quantum computation. We have observed 
similar shifts in other contexts which are discussed elsewhere [13]. 

Two important differences emerge in the modeling of external and internal er- 
rors. First, traditional boson-bath models[14] or spin-bath models[15], in which 
self-interactions among bath modes are omitted, are inappropriate for systems 
with important internal bath dynamics. While some master equations[16,17] 
include these effects, they are approximate and thus cannot be totally trusted. 
Second, the assumption of coordinate type couplings in boson-bath models re- 
sults in vanishing canonical averages for the bath coupling operator, and so 
no coherent shift is observed. Fortunately, the problem of internal errors in a 
small QC is amenable to exact propagation. We can thus avoid all the pitfalls 
of the approximate theories. 



The organization of this Letter as follows. In Section 2 we describe our isolated 
QC model and our exact numerical approach. In Section 3 we present our 
results for two error quantifiers, purity and fidelity, for a number of initial 
states and two different error generators. Section 4 discusses the observed 
effects in more detail. 



2 Isolated flawed QC model 

The total Hamiltonian of the isolated flawed QC reads 
H(t) =Hs(t) + HsB + H B 



in which Hs(t) is control Hamiltonian of the system (active part), Hb is the 
bath Hamiltonian (idle part), and Hsb is the interaction Hamiltonian. 
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Table 1 

The switching intervals and active Hamiltonians used to implement CNOT gate. 



Switching Intervals 



Active Hamiltonian 



[70 = 0,71=^/(2^)] 

[t 1 ,t 2 = t 1 + it/(2B x )} 

[t 2 ,T3= 72 +7T/(2£P)] 




[7-3,7-4 = 73 + ^/(2^)] 
[T 4 ,r 5 = T 4 + 7r/( 4 ^)] 

[r 5 ,r 6 = r 5 + y/2*/(2B x )] 



[t 6 ,t 7 = t 6 + tt/(2B z )} 
[t 7 ,t 8 = t 7 + it/(2B*)] 
[t$,t 9 = t$ + tt/{2B z )} 




The system Hamiltonian to control two active charge-qubits [12] is 



(2) 



i=i 



The CNOT protocol given in Ref. [12] can be implemented from Eq. (2) in 
nine steps with fast Hadamard gates, and this is summarized in Table I. Note 
that ignoring all possible imperfections in (2), and assuming that consecutive 
gates can be simultaneously switched on and off ( i.e. no free Hamiltonian 
evolution is allowed) we eliminate all likely systematic unitary errors in the 
system. (See Ref. [18] for a classification of errors). Hence, all unitary errors 
we observe will arise from interaction with the bath qubits. 

Defining two bath interaction operators T, x = K:^x an d ^ 2 = Yli!^ 2 ^z&l, 

the interaction Hamiltonian takes one of the two forms 



where a G {x, z}. Note that while the bit-flip errors by xx-type coupling occur 
naturally, we do not expect phase-errors by zz-type coupling to appear in this 
design. Nonetheless, we also explore phase-errors for pedagogical reasons to 
demonstrate how the decoherence and shift can vary for different system- 
environment interactions. 

Including one- and two-body static imperfections, the two-qubit Hamiltonian 
of the system (2) can be generalized for an N-qubit bath: 




(3) 
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-i iV+2 N+l N+2 

Hb = -^Y: (B*ai + B?*i) + E E ( 4 ) 

i=3 i=3 j=i+l 

Here, individual weights Bf are sampled randomly and uniformly, i.e. Bf € 
[Sq — 5/2, -Bq + 5/2]- -Bq is the average value of the distribution and 5 is 
a detuning parameter. The two-body residual interactions J* J and residual 
system-bath couplings X l a are also randomly and uniformly sampled within 
[—J, J], and [—A, A], respectively. Experimentally accessible control param- 
eters for the charge-qubits[12] are B a = 1.00 e, J x = 0.05 e and kT = 0.25 e, 
in units of e = 200 mK. The total gate time is then r 9 = 1.129 x 10~ 9 s. 
Since the idle and active qubits belong to the same QC, the idle qubits should 
only differ from active qubits by imperfections. Hence, the average of the 
distribution corresponds to Bq = B a , and we set 5 = 0.4 e. We consider 
values J x = 0.05, 0.25, 0.50, 1.00, 2.00, with units in e, to explore the on- 
set of chaos and determine whether the resulting bath-chaos is beneficial as 
predicted[6,7,9,10]. Note that the spectrum is Wignerian above J x = 0.15 e. 
(Exact diagonalization of the bath Hamiltonian is accomplished by a Lanczos 
algorithm[20] for iV = 10 qubits.) 

In order to avoid error propagation due to initial residual entanglement be- 
tween system and bath qubits, we assume an initial state of product form 

p(0) = p s (0) ® p B (0) (5) 



where ps(0) is the density of the active system qubits and ps(0) is the density 
of the idle bath qubits. In our simulations we consider eight different initial 
system states ps(0) = |V'o)(V'o| : f° ur standard basis states 

|Vo>e{|00>,|01),|10>,|ll)}, (6) 



and four Bell states 

|^ )G{(|00)±|ll))/v / 2,(|01)±|10))/v / 2}. (7) 



The assumption that the bath qubits are idle, and the fact that the residual 
two-body imperfections exist in the idle part, leads to dynamical thermaliza- 
tion of the bath qubits and hence the appropriate initial state is of canonical 
form. 

The reduced density of the system qubits is calculated as a partial trace over 
the bath degrees of freedom, 

ps(t) = Tr B [U(t)p(0)U\t)} (8) 
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where the propagator is 

t 

U(t) = texp [-(i/h) j H(t')dt'} (9) 



and T denotes the time ordering operator. The reduced density at time t can 
be expressed exactly as 

n « 9 -E n /kT 

ps(t) = E —ft- i^[|*„(f)K*n(t)l] (io) 

n=l V 



where the partition function is 

^eig 

Q = J2 e- E - /kT (11) 

n=l 

and H B \4>B) = E n \(p^) : and \^ n (t)) is a solution of the Schrodinger equation 
d\* n (t))/dt = -(i/h)H(t)\* n (t)) (12) 

with initial condition |\P n (0)) = \ipo) <8> Solutions of the Schrodinger 

equations were obtained using a variable stepsize Runge-Kutta code [21]. The 
number of initial states which are thermally populated is small at low tem- 
peratures and so only a few need be propagated. 



3 Results 



Two measures of errors were employed to quantify errors during the CNOT 
gate: average purity, and average fidelity. For each we compute two averages: 
over the four initial standard basis states, and over the four initial Bell states. 
The average purity, 

^ IV>0> 



measures proximity of the reduced density to a pure state. In the absence of 
system-bath coupling P(t) — 1 at all times for pure initial states. The average 
fidelity, 

m = lT,^s[p s m deal (t)] (M) 

IV'o) 
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measures how close the computed result is to the desired algorithm. In the 
absence of system-bath coupling, jF(t) = 1, and ideal time evolution is given 
by 

pf al (r 9 ) = t>cNOTPs(0)^ NO T- (15) 



In Fig. (la) we plot average purity vs time, and in Fig. (lb) average fidelity 
vs time for bit-flip errors for five different values of intra-bath couplings J x . 
Similar quantities are plotted in Fig. 2 for phase errors. In the Figures, the 
Bell states are distinguished from standard basis states by thick lines. Non- 
negligible deviations from perfect purity, greater than theoretically acceptable 
limit [23] of 0.99999, are observed in almost all cases, indicating that internal 
decoherence exists even for this relatively small number of bath qubits. How- 
ever, with increasing J x , a transition to chaos occurs, which in turn results 
in rapid suppression of non-unitary errors. The strongest intra-bath coupling 
J x = 2.00 e leads to nearly complete suppression of decoherence, which is 
true for both bit-flip and phase errors. Hence, deliberately induced bath chaos 
seems beneficial, and may even serve as error correction strategy when such 
strong interactions are practical. As expected, Bell states are more fragile to 
the effects of internal imperfections, and non-unitary bit-flip errors are of the 
same magnitude as phase-errors. 
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Fig. 2. (Color online) (a) Average purity P(t), and (b) fidelity J-{t) vs. time for 
phase errors. 



The fidelity plots show a very large and unexpected deviation from ideality 
even though the purity does not show this behavior. Hence, the large dis- 
crepancies between the purity and fidelity plots indicate that the environment 
primarily induced unitary errors via coherent shifting processes. Increasing 
J x results in an increase in fidelity for phase errors, but no-improvement is 
observed for bit-flip errors. The surprising magnitude of the observed unitary 
errors destroy practical QC operations, and methods need to be devised for 
incrementally correcting the dynamics so that multi-gate algorithms are vi- 
able. Note also that this error, here observed for a two-qubit subsystem, is of a 
collective nature and hence common to the entire active part which in general 
may consist of many qubits. While existing error correction strategies could 
be employed to remove this shift [19], the number of single-body operations 
would have to greatly increased for each qubit. Since the shift we observe is 
unitary and of a collective nature a much simpler collective error correction 
scheme may be possible. 



We performed our simulations on ten different realizations of the QC. The 
results we report here are typical. However, we encountered exceptions. Acci- 
dental near degeneracies in the low lying bath energies can occasionally occur 
with increasing J x , causing increased decoherence. Architectures involving xy- 
type two-qubit imperfections actually favor this behavior. In all cases, the 
large magnitude of the coherent shift is universal. 
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4 Discussions 

The observed reduction of decoherence with increasing J x can be explained by 
semiclassical arguments [9, 2 2]. Basically, the squares of the off-diagonal matrix 
elements of the active-idle coupling operator - in the basis of bath eigenstates 
- scale as fr^ -1 when the bath is chaotic[22]. For small H they therefore vanish 
in the thermodynamic limit. When this is combined with the fact that the 
diagonal elements vary slowly with energy in the low energy region of the 
spectrum, it guarantees that decoherence and dissipation will be less than 
would be found in an integrable bath where selection rules are operative. 

The coherent shift is a well known feature of the exact Nakajima-Zwanzig 
master equation[24], which is approximately of the form 

43s(t) = -(i/h)[H s (t) + SE,p s (t)] + Cpsit) (16) 



in the Markovian limit. (See also our detailed discussion of a similar shift in 
[13].) Here S = (d-^ + a^) is the system coupling operator and E = Tr[E a ps(0)] 
is the canonical average for the bath coupling operator. C includes non-unitary 
contributions to the dynamics. Clearly, the coherent shift originates from the 
dressing term SE in Eq. (16). In Fig. 3 we plot |E X | and |E 2 | as a function 
of J x . Increasing J x results in a small amount of increase in magnitude of 
|E X |, but a considerable amount of decrease in magnitude of |E 2 |, which is in 
very good agreement with the amount of decrease in fidelity for bit-flip errors 
in Fig. 1(b) and increase in fidelity for phase errors in Fig. 2(b). Chaos thus 
affects the shift indirectly through canonical averages which change with the 
bath eigenstates (i.e. E varies with J x through ](/>%))■ Since bit-flip and phase 
error generators do not commute, each has a different E even for the same 
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bath states, and thus their behaviors differ. It should also be noted that since 
increasing bath chaos results in improved purity, a corresponding improvement 
in fidelity is expected. However, unitary errors due to the shift are so strong 
that fidelity improvement is not noticeable. 

The fact that the coherent shifting is sensitive to intra-bath coupling strength 
is suggestive. First, it means that the shifting is at least partly a collective 
phenomenon. Second, while the shifting is clearly undesirable in a QC it may 
actually be useful in other situations. Note that Eq. (16) predicts that no 
shifting should be observed for a system coupled to a harmonic oscillator bath 
via a coordinate type coupling. In cases where shifting is observable it then 
gives a measure of anharmonicity and internal coupling of the bath degrees of 
freedom. Shifts observed in different types of optically active impurities could 
thus be employed to obtain detailed information about a solid. 



5 Conclusion 

In summary, the external decoherence time is not the only factor which should 
be taken into account in QC design and error correction strategies. Qubits in 
flawed QCs are subject to destructive effects of internal decoherence, dissipa- 
tion, and coherent shifting. Bath-chaos may serve to correct non-unitary errors 
but the unitary errors due to the coherent shifting remain. While existing er- 
ror correction codes may cure unitary errors, this will greatly increase the 
number of required one-qubit rotations. We hope that our preliminary results 
promote further study of architecture specific internal errors and strategies 
against them. 
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